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Spontaneous fission barriers, quadrupole inertia tensor, and zero-point quadrupole cor- 
relation energy are calculated for 252 ' 256 > 258 Fm in the framework of the self-consistent 
Skyrme-Hartree-Fock+BCS theory. Two ways of computing collective inertia are em- 
ployed: the Gaussian Overlap Approximation to the Generator Coordinate Method and 
cranking ansatz. The Skyrme results are compared with those of the Gogny-Hartree- 
Fock-Bogolyubov model. 

1. Introduction 

The microscopic description of the fission process is one of the most challenging 
issues in nuclear structure. It represents an extreme example of the tunneling of 
the many-body system. There exist many descriptions of spontaneous fission based 
on the adiabatic assumption. While the microscopic-macroscopic models still offer 
the best fit to the data, during recent years, a number of self-consistent approaches 
rooted in the Adiabatic Time-Dependent Hartree-Fock method (such as the collec- 
tive Schrodingcr equation with microscopic mass tensor) have been introduced. 

The spontaneous-fission half-life can be reliably estimated within the semi- 
classical WKB approximation. Here, the key quantity is the action integral 

along the trajectory s in the multidimensional space of collective parameters {q}. 
In Eq. (JTJ), V(q) represents the collective potential given by the HF+BCS or HFB 
energy, E is the ground-state energy and B(q) is a second rank tensor of a collective 
mass. 

The total energy of the system as a function of collective coordinates represents 
a zero-order approximation to the potential energy surface for shape- vibrations, 
V(q). However, before one can use V(q) in calculations with the collective Hamil- 
tonian, dynamic corrections have to be added. The reason is that the underlying 



9, 2008 0:46 WSPC/INSTRUCTION FILE parmasl3 



2 A. Baran, A. Staszczak, J. Dobaczewski, W. Nazarewicz 

states have a finite uncertainty in the collective deformation. As a consequence, the 
potential V(q) contains contributions from quantum fluctuations, and these con- 
tributions need to be subtracted first before adding the energies associated with 
the true physical zero-point fluctuations, Eq, see RefP The theoretical evaluation 
of these correction terms can be done in the framework of the cranking approxi- 
mation (CRA), the generator c oordinate method (GCM) P or the Gaussian overlap 
approximation to GCM (GOA)!^15lSEl 

Having minimized the action S in the collective space {q}, the spontaneous- 
fission half-life T s f can be calculated frorrP 

logT sf = -20.54 + log [1 + cxp (2S min )] - log (2E ). (2) 

It immediately follows from Eq. {1} that both the potential V{q) and the inertia 
tensor B(q) determine the fission dynamics in the multidimensional collective space; 
hence the value of T s f. 

In this work we discuss fission barriers, mass parameters, and zero-point correla- 
tion energy (ZPCE) calculated in the Hartree-Fock+BCS model with the commonly 
used Skyrme SkM* and SLy4 energy density functionals. The Skyrme re sults are 
compared to those within the Gogny D1S Hartree-Fock-Bogolyubov modelP^ An- 
other comparison is done between collective inertia obtained in the GOA and CRA 
approximations. 

2. The Model 

In order to make a comparison between various methods and models transparent, 
in the present calculations we only consider one collective variable q, an average 
axial quadrupolc moment of the mass distribution in the nucleus: 

g = Q 20 = (0|g 2O |0), (3) 

where 

A I A 

g 2 o = W nr^&'fc) = £( 2 *? - x ' - »?) • ( 4 ) 

i=i i=i 

To account for pairing correlations, we applied the BCS approximation to the 
density-dependent 5-pairing interaction (DDDI): 

u P air(n - r 2 ) = v ot S(ri - r 2 ) 

where t=p/n for protons/neutrons, p = p{r) is the isoscalar nucleonic density, and 
po=0.16 im~ 3 . The coupling constants «o n =842 MeVFm 3 and w Op =1020 MeVfm 3 
were fitted to the experimental pairing gaps in 252 Fm: A„ = 0.696 and A p = 
0.803 MeV. The pairing-active space consisted of the lowest Z (N) proton (neutron) 
single-particle states. We have also performed calculations with the seniority pairing 
force (with the constant strength parameters) given by 

G n = [24.70 — 0.108 (N — Z)] J A , 

G p = [14.76 + 0.241 (N - Z)\ /A. 1 ' 



1 - 



p(n) 



(5) 
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The calculations were performed using the code hfodd that 
allows for an arbitrary symmetry breaking. For the basis, we took the lowest 1140 
single-particle states of the deformed harmonic oscillator. This corresponds to 14 
oscillator shells at the spherical point. Based on the hfodd self-consistent wave- 
functions, the collective mass tensor components and ZPCE were computed. In prin- 
ciple, all the collective-tensor components A/i with multipolarities A = 2,3, ... ,9 
can be calculated. 

In this study, only specific static fission paths are considered; they have been ob- 
tained in the calculations presented in RefP^ Specifically, for the fermium isotopes 
considered, one predicts both reflection-symmetric (s) and reflection-asymmetric 
(a) fission valleys. There are two kinds of symmetric paths predicted; namely, the 
valley that corresponds to elongated fission fragments (E), and that with more com- 
pact (C) fragments that resemble spherical 132 Sn clusters when approaching 264 Fm. 
(The shorthand notation sEF, used in the following, means symmetric, elongated- 
fragments fission path; similarly aEF and sCF.) 



3. Inertia Tensor and Zero-Point Correlation Energy 

In order to calculate the inertia tensor one usually applies the GOA or uses the 
cranking model formalism. In both cases, the main task is the calculation of M^ k \ 

mW (n) v (m\^)(^\Q}\^ (7) 

vv' 

where qt = (Qi) are the collective coordinates (i stands for Xfi), \vv') = a£,aj,, |0) 
is a two-quasiparticle state, and E v is the quasiparticle energy. 

Both i n th e GOA and cranking, the inertia tensor can be given by compact 
expressions^^ 

B GOA = E (2)[ S (l)]-l E (2) ) (8) 

B CRA = S (3) ^ (Q) 

where matrices S^- 1 read 

sI^V'-'mWm' 1 )- 1 , (10) 

4 v ' 

In the nuclear case two kinds of fermions are present. The total covariant inverse 
inertia for a composite system is given as a sum of proton and neutron covariant 
inertia tensorsP This leads to the final expression 

a B n B 

IpBn + ltB p 

where 7 is the total metric tensor, which is a sum of proton (7^) and neutron (jn) 
contributions. 
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Fig. 1. The total (t) GOA quadrupolc collective inertia and its proton (p) and neutron (n) com- 
ponents calculated for 258 Fm in the SkM* model with seniority pairing along the asymmetric 
elongated static path to fission. 

We conclude this section by recalling that the zero-point correlation energy in 
the GOA can also be expressed through quantities (fT0|) !^ i.e., 

E a = ^TrCE^- 1 ^). (12) 

4. Results 

We begin our discussion from Fig.[TJ which shows the quadrupole GOA inertia B 2 o 
© for 258 Fm calculated along the aEF path in the SkM* -(-seniority model. The 
corresponding proton and neutron contributions are also depicted; they are related 
to the total inertia through Eqs. (flT)|) and (fTTj) . The behavior of the total inertia can 
easily be traced back to the proton and neutron shell effects that are responsible 
for characteristic fluctuations in the collective mass. 

Figure [2] compares the values obtained in SLy4+DDDI and D1S models for 
252 Fm along the sCF static fission path. While the SLy4 values usually exceed those 
obtained with the Gogny force by more than a factor of two. the general patterns 
of -B2o(<?2o) ar c fairly similar. The same conclusion holds for comparison between 
GOA and cranking inertia, with the cranking value being larger. A very interesting 
feature of the self-consistent inertia parameters is their regular behavior at large 
elongations. This has not been seen in earlier calculations using phcnomcnological 



9, 2008 



0:46 WSPC/INSTRUCTION FILE parmasl3 



Collective inertia and fission barriers within the SHF theory 5 



CRA-D1S 




Fig. 2. Comparison of B20 in GOA and CRA for 252 Fm calculated along the symmetric compact 
path to fission in Skyrme-HF (SLy4+DDDI) and Gogny-HF (D1S) models. 

potentials (see, e.g., Ref.pl) where the oscillating behavior of B20 persisted in the 
whole deformation range. 

In order to calculate the collective potential V(q) and fission barriers, the zero- 
point correlation energy has to be evaluated first. The ZPCE (fT2"j) for 252 Fm cal- 
culated in SLy4+DDDI and D1S is shown in Fig. [3l One sees a qualitative and 
quantitative agreement between both variants of calculation. The largest correc- 
tion is predicted around the first minimum (Q20~25 b). As seen in Fig.[4l the inner 
and outer fission barriers in 252 Fm arc significantly modified by the ZPCE. Namely, 
the barriers calculated with the corrected collective potential are ~lMeV higher in 
the whole range of Q20- 

Finally, Figs. [5] and [6] display the collective potentials and inertia parameters (in 
GOA and CRA) in 256 Fm and 258 Fm, respectively, calculated along the symmetric 
(sEF and sCF) and asymmetric (aEF) fission pathsP^l These results confirm the 
previous observations; namely, the deformation patterns of collective inertia in CRA 
and GOA are very similar, with the cranking values being appreciably higher. 

5. Summary 

In this study we performed pilot calculations of collective inertia, zero-point 
quadrupole correlation energy, and fission barriers in selected Fermium isotopes. 
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Fig. 3. The zero-point correlation energy i ] 1 21) in GOA along the aEF fission path in 252 Fm 
calculated in SLy4-DDDI (+) and D1S (X) models. 



> 

CD 



o 

>s 

O) 
CD 

LU 

-t— • 

o 




50 100 150 

Quadrupole Moment Q 20 (b) 



200 



Fig. 4. The collective potential in 252 Fm in SLy4-DDDI along the aEF fission path without (E tot ; 
thin line) and with ( E*°l - E ; thick line) the ZPCE included. 
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Fig. 5. The collective potential (left-hand-side scale) and corresponding GOA and CRA quadrupole 
inertia (right-hand-side scale) along three fission paths in 256 Fm using the SkM*+seniority model. 



The hfodd code employed allows for an arbitrary symmetry breaking; this feature 
is of crucial importance when discussing spontaneous fission where the reflection- 
asymmetric and triaxial shapes can play a role. The main conclusion of this work 
can be summarized as follows: 

• The zero-point correlation energy is important to include as it can signif- 
icantly modify the action by lowering the ground state by approximately 
1 MeV in the nuclei considered. 

• The collective quadrupole inertia calculated in the Skyrmc SLy4 and SkM* 
models have very similar deformation pattern to those in th e Gq gny D1S 
model. Our HF-Skyrme values arc very close to those of Refs!-^ 

• The collective inertia obtained in the cranking approximation are about 
twice as large as the GOA results. In all the cases considered, the collec- 
tive inertia smoothly decrease at large deformations; they do not exhibit 
oscillations seen in previous microscopic- macroscopic calculations. 

• As discussed earlier, both V(q) and B(q) contribute to the action integral 
S (IH). In general, the fission barriers in the Gogny model arc higher as 
compared to those in SkM* and the opposite holds for the collective inertia. 



9, 2008 0:46 WSPC/INSTRUCTION FILE parmas!3 



8 A. Baran, A. Staszczak, J. Dobaczewski, W. Nazarewicz 




50 100 150 200 250 300 350 400 450 
Quadrupole Moment Q [b] 



Fig. 6. Similar to Fig. [5] except for 258 Fm. 



However, both those effects partly cancel in the product B(V — E) that 
determines the action and the spontaneous fission lifetime. 
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